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Probing intergalactic radiation fields during cosmic 
reionization through gamma-ray absorption ^ 
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ABSTRACT 

We discuss expectations for the absorption of high-energy gamma-rays by 77 pair 
production with intergalactic radiation fields (IRFs) at very high redshifts (z ^ 5 — 20), 
and the prospects thereof for probing the cosmic reionization era. For the evolving IRF, 
a semi-analytical model incorporating both Population II and Population III stars 
is employed, which is consistent with a wide variety of existing high-z observations 
including QSO spectral measurements, WMAP Thomson depth constraints, near-IR 
source count limits, etc. We find that the UV IRF below the Lyman edge energy 
with intensities in the range of a few times 10~^^ erg cm~^s~^Hz~^sr~^ can cause 
appreciable attenuation above '^12 GeV at z 5, down to 6 — 8 GeV at z ^ 
8 — 10. This may be observable in the spectra of blazars or gamma-ray bursts by the 
Fermi Gamma-ray Space Telescope, or next generation facilities such as the Cherenkov 
Telescope Array, Advanced Gamma-ray Imaging System or 5@5, providing invaluable 
insight into early star formation and cosmic reionization. 

Key words: galaxies: high-redshift - intergalactic medium - cosmology: theory - 
gamma-rays: bursts - galaxies: active 



1 INTRODUCTION 

Some time after the epoch of cosmic recombination at red- 
shift z ~ 1100, the bulk of the intergalactic gas in the 
universe must have been somehow reionized by 2 ~ 6, as 
indicated observationally from the spectra of high-2 QSOs 
and the polarization of the cosmic microwave background 
(CMB). However, the sources, history and nature of this cos- 
mic reionization process are still largely unknown, as most 
of this redshift range has yet to be explored through di- 
rect observations. Because the first stars and galaxies in the 
universe must have formed during this period, the primary 
suspect is photoionization by UV radiation from such ob- 
jects, potentially involving metal-free. Population (Pop) III 
stars. Alternative possibilities include mini-quasars, super- 
nova remnants and dark matter decay. Besides providing us 



* Numerical data of the model results will be available at 
http: //www-tap. scphys.kyoto-u.ac.jp/~inoue/hizabs/ 
t E-mail: inoue@tap.scphys.kyoto-u.ac.jp 



with clues to such processes in the early universe, cosmic 
reionization also profoundly affects the ensuing formation of 
stars and galaxies, so elucidating this era is one of the most 
pressing issues in cosmology today (see Barkana & Loeb 
2001; Ciardi & Ferrara 2005; Fan et al. 2006, Choudhury 
2009 for reviews). 



In the majority of scenarios for reionization of hydro- 
gen in the intergalactic medium (IGM), the main protago- 
nists are UV photons with energies above the Lyman edge 
(e > CLE = 13.6 eV). Although those with lower energies 
do not contribute to photoionization, they are also crucial 
since i) they give indications as to the strength and nature of 
the ionizing radiation, ii) those in the Lyman- Werner band 
(e = 11.2 — 13.6 eV) can photodissociate H2 molecules and 
suppress early star formation (e.g. Ciardi & Ferrara 2005), 
and iii) Lya photons (e — 10.2 eV) can strongly affect the 
HI spin temperature and the associated cosmological 21 cm 
signatures (e.g. Furlanetto et al. 2006). Thus, having some 
observational means to probe the evolution of UV intergalac- 
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tic radiation fields (IRFs) ^ in the cosmic reionization era 
would be of paramount importance, complementing exist- 
ing observations that probe the neutral or ionized gaseous 
components of the IGM. However, direct detection of this 
diffuse emission from very high z is extremely difficult if not 
impossible. ^ 

An indirect but powerful means of probing diffuse radia- 
tion fields is through photon-photon (77) absorption of high- 
energy gamma-rays (e.g. Gould & Schreder 1967; Stecker et 
al. 1992). Gamma-rays with energy E emitted from extra- 
galactic sources will be absorbed during intergalactic prop- 
agation by interacting with photons of the diffuse radia- 
tion field with energy e to produce electron-positron pairs 
(7 + 7 ^ e"*" -|- e~), as long as there is sufficient opacity for 
energies satisfying the threshold condition Ee{l — cosO) > 
2rnj.c* , where 9 is the incidence angle of the two photons. 
The observed spectra of the gamma-ray sources should then 
exhibit corresponding attenuation features, from which one 
can effectively infer or limit the properties of the diffuse ra- 
diation. This method has been utilized in recent TeV obser- 
vations of blazars by ground-based Cherenkov telescopes to 
set important constraints on the cxtragalactic background 
light in the near infrared to optical bands at relatively low 
z (Aharonian et al. 2006; Albert et al. 2008). 

As first discussed by Oh (2001; see also Rhoads 2001), 
UV IRFs with sufficient intensities to cause IGM reioniza- 
tion are also likely to induce significant 77 absorption in 
gamma-ray sources at ^ > 6 at observed energies in the 
range of a few to tens of GeV. However, these estimates i) 
were made before WMAP observations indicating an early 
start of reionization and were limited to 2 < 10, and ii) 
did not include the possibility of metal-free Pop III stars, 
which may have been active during the first epochs of star 
formation and are more prodigious UV emitters compared 
to normal stars ^. The recent launch of the Fermi satellite 
II with the Large Area Telescope (LAT) operating in the 
~ 0.1 — 100 GeV domain motivates us to reevaluate the 77 
absorption opacity at very high z, incorporating more recent 
observational and theoretical developments concerning the 
cosmic reionization era. 

For this purpose, we employ updated versions of the 
semi-analytical models of Choudhury & Ferrara (2005; 
2006), which self-consistently describe inhomogeneous reion- 
ization of the IGM, accounting for both Pop II and Pop III 
stars and their radiative and chemical feedback effects. With 
only a few free parameters, they axe able to fit a wide var 
riety of high-0 observational data. Using the evolving IRFs 



t Although often referred to as "extragalactic background light" 

for lower 2, here we avoid the term "background", since IRFs can 

be highly inhomogeneous in the reionization era, especially for 

e > CLE, even though it turns out to be more or less uniform for 

the spectral regime relevant to 77 absorption; see below. 

§ Earlier indications of a large contribution from Pop III stars 

to the local near-IR background are now disfavored from J-band 

source count limits (Salvaterra & Ferrara 2006). 

^ Note that 77 absorption measurements in \ow-z blazars have 

set strong constraints against a large contribution from Pop III 

stars to the local near-IR background (Aharonian et al. 2006; see 

also Raue et al. 2009). 

II http: / / fermi.gsfc.nasa.gov 



as predicted by these models, the 77 opacity is evaluated 
for the redshift range 2 = 5 — 20. We also briefiy assess the 
detectability of the resultant absorption features in high-« 
sources such as blaaars or gammar-ray bursts (GRBs) with 
current and future gammar-ray facilities, and the consequent 
implications. 



2 INTERGALACTIC RADIATION FIELD 
MODEL 

The salient features of our semi-analytical models are as fol- 
lows (sec Choudhury & Ferrara 2005; 2006; 2007; Choud- 
hury et al. 2008; Choudhury 2009 for more details): (1) 
Adopting a lognormal distribution of IGM inhomogeneities 
(Miralda-Escude et al. 2000) , the ionization and thermal his- 
tories of the neutral, HII and Hell phases of the IGM are 
tracked siinulaiieously and sclf-consistcntly. (2) The fornia- 
tion and evolution of dark matter halos are described by a 
Press-Schechter-based approach. (3) Three types of radia- 
tion sources are considered: a) metal-free Pop HI stars with 
a Salpeter initial mass function (IMF) in the mass range 
1 — IOOM0, with spectra according to Schaerer (2002) and 
including nebular and Lya emission lines (see Salvaterra & 
Ferrara 2003); b) low-metalhcity {Z = 0.02^©) Pop II stars 
with spectra according to Bruzual & Chariot (2003), other- 
wise being the same as Pop III; and c) QSOs with power- 
law spectra and emissivity based on the observed luminosity 
function at 2 < 6, considering only those above the break 
luminosity (Choudhury et al. 2008). (4) Pop II and Pop HI 
stars each form from gas in virializcd halos with efficiencies 
e*,ii and e*,iii, respectively, and the corresponding escape 
fractions of ionizing photons from the host halos are pa- 
rameterized by /esc,ii and /esciii- Included self-consistently 
are the consequent effects of radiative feedback that sup- 
presses star formation in sufficiently small halos, as well as 
a "genetic" merger-tree-b£ised treatment of chemical feed- 
back that induces the transition from Pop HI to Pop II star 
formation (Schneider et al. 2006). 

The free parameters of the model are e*,ii, e,,iii, Tjesc, 
which fixes both /esc. 11 and /osc.iii, and Ao, related to the 
mean free path of ionizing photons due to HI in high-density 
regions. These are ascertained so as to simultaneously re- 
produce a largo set of high-2 observational data: i) redshift 
evolution of Lyman-limit absorption systems; ii) effective 
optical depths of the IGM for Lya and Ly/3 from QSO spec- 
tra; iii) electron scattering optical depth re from WMAP 

3rd year results (Spergel et al. 2007); iv) temperature of 
the mean IGM; v) cosmic star formation history; and vi) 
limits on J-band source counts from NICMOS HUDF. In 
the fiducial, best-fit model, H reionization begins rapidly 

at 2 ~ 15, initially driven by Pop III stars, and is 90 % 
complete by 2 ^ 8. Thereafter it is slowed down by feedback 



** The adopted cosmological parameters are h = 0.73, Clm = 
0.24, Ha = 0.76, Cl^h'^ = 0.022, as = 0.74, Us = 0.95 and 
dus/dlnk = (Spergel et al. 2007). 

tt The fiducial model gives Te = 0.07, consistent with the 5th 
year results as well (Dunkley et al. 2009). 

tt The relevant parameters are £« = 0.1, £« = 0.02, /esc, 11 = 
0.0578 and /esc,iii = 0.54. 
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Figure 1. Volume- averaged intensity of the intergalactic radia- 
tion field J(e) vs. energy e (or wavelength A) at redshifts z as 
labelled for the fiducial model. 



effects and taken over by Pop II stars at z ~ 7, finally reach- 
ing completion by z; ~ 6 (see Fig.2 of Choudhury 2009). The 
cosmic star formation rate is always dominated by Pop II 
stars and is at the level of 0.05 - 0.08 M© yr"^Mpc"^ for 
z ~ 6-8 (Fig.2(b) of Choudhury 2009), in hue with that de- 
duced from observed GRB rates (e.g. Salvaterra et al. 2008, 
Kistler et al. 2009). 

Shown in Fig.l is the volume-averaged intensity of the 
IRF as calculated from this model, which declines monoton- 
ically with z following the evolution of the star formation 
rate (SFR). We caution that at z > 6 before intergalactic 
HII regions have completely overlapped, the IRF is expected 
to be inhomogeneous and fluctuating along different lines 
of sight, particularly strongly for e > ele. However, it is 
also evident that ionizing photons are strongly absorbed by 
the neutral IGM and the mean IRF spectrum cuts off very 
sharply above ele, so that this portion has negligible effects 
on the 77 opacity (Madau & Phinney 1996, Oh 2001, Rhoads 
2001). On the other hand, UV radiation with e < ele have 
much longer mean free paths in the IGM, and the notion of a 
nearly uniform and isotropic background may still be appro- 
priate for this regime, which is also the most relevant for 77 
absorption. Of particular note is the band e — 10.2 — 13.6 eV 
where the spectrum dips somewhat due to blanketing by the 
Lyman series lines, but which should nevertheless be very 
important for the 77 opacity. 

Being optimized for the cosmic reionization era, the 
main shortcoming of the present model is that it does not 
account for Pop I stars or dust that can become important 
at lower z. Our IRF calculations are available only for z > 4, 
and may be somewhat less reliable near z ~ 4 as the com- 
parison with observations has not been as thorough as for 
z > 5. A more complete model describing the evolution of 
the IRF at all redshifts awaits future studies. 




Figure 2. Local 77 optical depth Tio^alC^') vs. rest- frame gamma- 
ray energy -Ercst at redshifts z as labelled for the fiducial model. 
The contribution from the CMB is not shown. 



3 GAMMA-RAY ABSORPTION OPACITY 

We first estimate the "local 77 optical depth" at each 
z by the optical depth across a Hubble radius Ih{z) = 

POO 

Tioz^i{z,E) = Ih{z) I den{e,z) 

^\ j "^/^(^ ^ ^)°"^^(-^''^'/^)' 

where n(e, z) is the IRF photon number density per energy 
interval, /i = cos0, eth = 2mgC*/_E(l — /i) is the thresh- 
old energy, and a^^{E,e, ^) is the 77 pair production cross 
section. For given E, a-y^ rises sharply from e — eth, peaks 
at e = 2eth, and then falls off as . Thus riocai roughly 
mirrors the IRF spectrum at each z, although its detailed 
features are smeared out. Displayed in Fig.2 in terms of the 
rest-frame gamma-ray energy E^cst, we see that the opac- 
ity may be significant out to z ~ 10 for -Brest ~ 10^ — 10* 
GeV. Note the steep drop in riocai at -Brest < Sle ~ 
(meC^)^/eLE — 18 GeV, corresponding to the sharp cutoff 
in the IRF spectrum above the Lyman edge. As pointed out 
by Oh (2001; see also Rhoads 2001), this is crucial in that it 
allows appreciable contributions to the total 77 opacity from 
higher z even when the IRF intensity is relatively weaker, 
and which should be uncontaminated from absorption at 
lower z. However, we also see that due to the declining IRF 
intensity together with the reduced path length, the opacity 
from z ^ 10 is likely to be quite small. 

This can be seen more explicitly in Fig. 3 where we show 
the integrated 77 optical depths for different source redshifts 



§§ The contribution from the CMB (e.g. Stecker et al. 2006) also 
becomes important at -Ercst > 2 TeV, but is irrelevant for our 
results below and not plotted in Fig.2. 
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Figure 3. Integrated 77 optical depth t{E) vs. observed gamma- 
ray energy E for source redshifts z as labelled, for the fiducial 
model (solid) and K04's high stellar UV model (dashed). 
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where dl/dz' = Ih{z') / {l-\-z') and E is the observed gamma- 
ray energy at z = 0. As mentioned above, the lower limit of 
^-integration that can be taken in our model is z-mi-n ~ 4; for 
additional absorption from the range z = — 4, we can only 
consult other models at the moment (e.g. Kneiske et al. 2004, 
hereafter K04; Stecker et al. 2006; Razzaque et al. 2009; 
Gilmore et al. 2009). Overlayed here for comparison is K04's 
"high stellar UV model" , which gives their best description 
of QSO proximity effect measurements at 2: ~ 2 — 4. 

As we could infer from Fig. 2, our model predicts appre- 
ciable opacity at observed energies < 12 GeV for sources 
at 2; > 5, with notable differences out to z ~ 8. However, 
the relative effects of further absorption from 2 J> 8 may be 
practically indiscernible. Nevertheless, the spectral attenu- 
ation feature itself should be observable in high-z gamma- 
ray sources by current or future gamma-ray facilities, and 
possibly distinguishable in the range z ~ 5 — 8 for suffi- 
ciently bright objects (§4). Owing to the drop in 77 opacity 
at Brest < 18 GeV (Fig. 2), the differences in absorption in 
this z range are caused in-situ by the evolution of UV IRFs 
just below the Lyman edge energy, including the crucial Lya 
and Lyman- Werner bands. We also recall that in this model. 
Pop 111 stars continued to be significant contributors to the 
UV IRF down to z ~ 7, where they are comparable with Pop 
11 stars for ionizing photons. Measurements of these effects 
would thus provide an important check of current models of 
cosmic reionization in its latter stages, as well as a unique 
and invaluable probe of evolving UV IRFs in the sub- Lyman 
edge regime during the era of early star formation (§4). 

In Figs. 4 and 5, respectively, we plot the spectral atten- 
uation factor exp[— r(i?)] and the observed energy E{t = 1) 
where the optical depth is unity. Here the fiducial results are 



compared with those of an alternative model that does 
not include Pop HI stars, in which reionization is driven 
only by Pop II stars and occurs relatively late at z ~ 6 
(similar to the late reionization model of Gallerani et al. 
2008). The fact that Pop II stars are less efficient sources 
of ionizing photons compared to Pop 111 stars mandates a 
larger SFR, more intense IRF for e < ele and hence larger 
77 opacity. However, since the SFR at 2 < 6 is observation- 
ally constrained, notable differences appear only at 2 ^ 8, 
which should be challenging to distinguish in practice. Thus 
gamma-ray absorption may not be a sensitive probe of the 
reionization history itself. We have also investigated various 
other models, e.g. those with more realistic prescriptions 
for radiative feedback that fit the current high-2 observa- 
tions nearly equally well, and found that they generally do 
not lead to large differences. Conversely, being constrained 
by existing data, our predictions may be considered reason- 
ably robust, at least within the framework of our model. 
Nevertheless, we caution that relaxing some of the present 
assumptions, e.g. regarding the stellar IMF or the QSO con- 
tribution, may yet allow a wider range of possibilities. Note 
that although some other recent models (e.g. Razzaque et 
al. 2009, Gilmore et al. 2009) predict somewhat less absorp- 
tion at 2 ~ 5 — 6, they are not directly comparable with ours 
as their focus is on the 2 < 6 universe (e.g. Gilmore et al. 
2009 do not attempt to fit the Lya effective optical depths 
at 2 ^ 5.5 as we do). 



4 DISCUSSION 

The fact that 77 absorption is sensitive to photons with 
energies below the Lyman limit rather than the ionizing ra- 
diation (§2) actually points to a unique probe of the cosmic 
reionization epoch that complements measurements of QSO 
Gunn-Peterson troughs or CMB polarization anisotropies, 
which probe the neutral and ionized components of the IGM, 
respectively. On the one hand, observationally deducing the 
global UV emissivity and hence the cosmic star formation 
rate from the latter two is problematic due to uncertainties 
in the inhomogeneity of the IGM (clumping factor) and the 
escape fraction of ionizing photons from the host galaxies 
(Madau et al. 1999, Wyithe et al. 2009). On the other, the 
direct census of the high-2 UV luminosity density from deep, 
near-lR surveys are affected by the uncertain integrated con- 
tribution of faint galaxies below the telescope detection limit 
(e.g. Bouwens et al. 2007). Observing 77 absorption in high- 
2 sources may allow more robust measurements of the evo- 
lution of the cosmic UV emissivity, and in combination with 
other data, possibly the determination of the escape frac- 
tion and/or the IGM clumping factor as well. These infer- 
ences are general and independent of any particular model 
for reionization, but will be investigated in more quantita- 
tive detail in the near future. Likewise, the implications for 
constraining Lya or H2-dissociating radiation from 77 ab- 
sorption will be discussed in future work. 

We now briefly address whether the effects discussed 



The relevant parameters are etji = 0.1 and /esc, II = 0.0928. 
The model gives Te = 0.06, marginally consistent with the 5th 
year WMAP constraints. 
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Figure 4. Spectral attenuation factor cxp[~T{E)] vs. observed 
gamma-ray energy E for source redshifts z as labelled. Upper 
panel: fiducial model (solid) and K04's high stellar UV model 
(dashed). Lower panel: model without Pop III stars. 
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Figure 5. Observed gamma-ray energy E where r = 1 vs. red- 
shift 2, for the fiducial model (solid), model without Pop III stars 
(dotted) and K04's high stellar UV model (dashed). 



above are observable in real sources witli current or future 
gamma-ray instruments. For blazars, tlie most prominent 
and numerous extragalactic sources of GeV gamma-rays, tiie 
highest redshift confirmed so far is z ~ 3 (Hartman et al. 
1999; Abdo et al. 2009b). However, objects similar to the 
most powerful known blazars such as 3C454.3 with apparent 
luminosities L ~ 10*^ erg s^^ should be detectable by Fermi 
out to 2 ~ 8 — 10 if they exist at such redshifts (e.g. Romani 
et al. 2004). According to the latest blazar evolution models 
(Inoue & Totani 2009), it may be plausible for Fermi to 
detect some blazars above z ^ 6 during its survey period. 



for which deep, pointed observations may indeed reveal the 
IRF absorption features described above. 

GRBs are also promising as they are known to occur 
at 2 > 6 (Kawai et al. 2005; Greiner et al. 2009), at least 
up to 2 8.2 (Tanvir et al. 2009; Salvaterra et al. 2009), 
and perhaps out to the very first epochs of star formation 
in the universe (e.g. Bromm & Loeb 2006; Salvaterra et 
al. 2008). Although the spectral properties of GRBs in the 
GeV domain are still rather uncertain, previous detections 
by CGRO/EGRET (Hurley et al. 1994) and the recent de- 
tection of GRB 080916C at 2 = 4.35 by Fermi (Abdo et al. 

2009a) III! demonstrate that at least some GRBs have lu- 
minous emission extending to >10 GeV, which can also be 
expected theoretically (e.g. Zhang & Meszaros 2001; Asano 
et al. 2009). A burst similar to GRB 080916C may stiU be 
detectable at several GeV by Fermi/LAT at 2 < 7, and 
even out to higher 2 if the spectrum was somewhat harder. 
Even better for this purpose would be proposed ground- 
based telescopes with much larger effective area and multi- 
GeV energy threshold, such as the Cherenkov Telescope 
Array (CTA)* * *, Advanced Gamma-ray Imaging System 
(AGIS)* * * or the 5@5 array (Aharonian et al. 2001). To- 
gether with measurements of Lya damping wings (McQuinn 
et al. 2009) and possibly radio dispersion (loka 2003; In- 
oue 2004), future, broadband observations of very high-2 
GRBs should open new windows onto the cosmic reioniza- 
tion epoch. 

Even if IRF-induced spectral features are detected, a 
generic problem for 77 absorption studies is distinguishing 
them from spectral cutoffs intrinsic to the source. In this 
regard, spectral variability should offer an important clue. 
Both blazars and GRBs are highly variable gamma-ray emit- 
ters, and in general, changes in physical conditions of the 
source that cause variations in fiux should also be accompa- 
nied by variations of the intrinsic cutoff energy, whether it 
is due to injection of freshly accelerated particles, changes 
in the magnetic fields, internal radiation fields, bulk fiow ve- 
locity, etc. In contrast, cutoffs of IRF origin should be sta- 
ble in time and independent of the variability state of each 
object. Acquisition of time-resolved spectra should thus al- 
low the deconvolution of the two effects. Another indication 
should come from statistical studies of a sufficient sample of 
measurements. IRF-related cutoffs should occur at similar 
energies for sources at similar 2, and also exhibit a system- 
atic evolution toward lower energies for higher 2, whereas 
there is no strong reason to expect such trends for intrinsic 
cutoffs. Both the above strategies motivate the construction 
of future, high-sensitivity multi-GeV facilities such as CTA, 
AGIS and 5@5, which should be powerful tools to probe the 
evolution of UV IRFs in the cosmic reionization era through 
77 absorption in very high-2 sources. 



Ill For z = 4.35 and E = 13.2 GeV, our fiducial model gives 
Tjj ~ 0.4, consistent with the actual detection of a photon at 
this energy from GRB 080916C. 

* * * http://www.cta-observatory.org 

* * * http://www.agis-observatory.org 
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